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Abstract 



Given a sequence /i(xi), f2{xi, X2), ■ ■ ■ , fk{xi, . . . , Xk) of Boolean functions, each of 
which fi takes the value 1 in a single point of the form x?, X2, ■ ■ ■ , x^, i = 1,2, . . . ,k. A 
length of all x^ is n, A = 2". It is shown how to find x^ {k > 2) using ^^^^ simultaneous 



evaluations of functions of the form fi, /j+i with an error probability of order k/\nV which 
is \/2 times as fast as by the k sequential applications of Grover algorithm for the quantum 
search. Evolutions of amplitudes in parallel quantum computations are approximated by 
systems of linear differential equations. Some advantage of simultaneous evaluations of 
all /i, . . . fk are discussed. 

1 Introduction 

1.1 Structure of the work 

After background and setting of a problem I present the short introduction to abstract quantum 
computations in the section 3. In the section 4 linear differential equations are applied to a 
tight analysis of the famous Grover algorithm of the fast quantum search. 

The section 5 is the key. Here a parallel quantum algorithm for repeated search is defined 
and studied by means of differential equations. In the section 6 we briefly run through a parallel 
algorithm for iterated search with simultaneous queries of all oracles. The abstract includes 
the other substantiation of the algorithm. 



1.2 Background 

The most promising quantum mechanical application to the algorithm theory is associated with 
the fundamental algorithm of exhaustive search or finding a solution of equation f{x) = 1 for a 
Boolean function /. In 1996 Lov Grover in the work ||Grl|| showed how quantum computer can 



solve this equation for the case of unique solution in a time O(v^) where N is the number of 
all possible values for x, whereas every probabilistic classical algorithm requires a time 0{N). 
At about the same time in the work ||BBBV|| it was shown that there are not substantially 



faster algorithm for this problem. Later a tight estimation for the time of Grover's algorithm 



*Department of applied mathematics, Moscow state university of technology "Stankin", Vadkovsky per. 3a, 
101472, Moscow, Russia, e-mail: y@ oz.msk.ru 



1 



as T^^^ with the probabihty of error about was estabhshed in the work [ BBHT |. A further 
development of the fast quantum search can be found in the works | |Dfd[ | , [ [b'G] ] , | |Jo| ], [ |CGW| ], 
IBBBGDLj , pMHj , 0. 
The earher patterns of quantum speedup was constructed by P. Shor (look at the work 



1 5^). There are the algorithms finding a factorization of an integer and a discrete logarithm. 
Algorithms of such a sort are presented in a lot of works (look for example at |DJ|, ||Gr2|| , ||Kil|] , 



t|, lTg,etc.). 

The classical computations admitting quantum speedup are rare exclusions from all classical 
computations in the following sense. Denote all words of a length n in the alphabet {0, 1} by 
{0, 1}". We can represent a general form of classical computation as T iterated applications of 
some oracle g : {0, 1}" — > {0, 1}": 



X — > g{x) — > g{g{x)) 



g{g{...{g{x)))). 



(1) 



In the work ||Ozl|| it is shown that if T = 0{N^), e > 0, then for the bulk of all g such 
computation has not any quantum speedup. Similar results for search problems were obtained 
in the works ||BBB V|| , [|BBHT|| , [Pz2|| , [^al . A lower bound as 0{N) for a time of quantum 
computations of functions with functional argument F : {/} — > {0, 1} was found in the work 



|BBCMW| . Here / are Boolean functions on domain of cardinality A^. At the same time using 
a memory of 0{N) qubits it is possible to compute such functions in a time A^/2 (look at [[vD|| ) . 
This brings up the question: what a general type of classical computations of the form (1) 
admits a quantum speedup beyond any possible speedup of g7 As follows from the work [ pzlf 
for the bulk of functions g this speedup can result only from a parallel application of g. 
About other aspects of quantum evolutions look also at | [ML| | , [ [Ho| | , | |Fe|| . 



2 Setting of the problem 

Consider the following situation. We want to gather a mosaic from scattered stones in a 
rectangular list with the corresponding picture. Each stone is of a unique form. We can gather 
this mosaic layer by layer and use a simple search among stones still scattered to fill any layer 
basing on the previous one. Then we in fact fulfill an iterated search classically, because to find 
the stones for the following layer we must already have the previous one filled. 

We formalize this as the special type of iterated algorithms: iterated search (IS). Suppose we 
have a sequence 5*1, 5*2, . . . , S'^ of similar search problems where Si is to find a unique solution 
x° of equation fi{xi) = 1 where a Boolean function fi is accessible iff we know all x^, j < i, 
\xi\ = n, N = 2", k N, \x\ denotes a length of word x. The aim is to discover k > 2. 
In view of the result of the work | |BBH'1' | cited above sequential applications of Grover's search 



for x^,X2, ■ ■ ■ give an answer in the time with error probability about k/N. To do 

this we must have all oracles ft, i = 1,2, . . . , k, where the dependence fi of all Xj, j < i can 
be included to fi. So we can assume that fi has the form fi{xi,X2, . . . ,Xi) and each equality 
fi{xi, . . . ,Xi) = 1 has the unique solution x?, Xg, . . . , i = 1,2, . . . , k. Considering all oracles 
fi as physical devices which can not be cloned we assume that they are in our disposal at 
the same time, so we can apply them simultaneously. In such application advantage is taken 
of interference between the results of their actions. This results in a speedup of computation 
comparatively with the sequential mode. Why this speedup can arise? It arises because of a 
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leak of amplitude in each step of sequential search. A leak of amplitude issues from that an 
amplitude of x° in search number i increases step by step in course of Grover's search, after few 
first / steps it becomes approximately when amplitudes of others Xi ^ decrease. This 
prevailing of the amplitude corresponding to (a leak of amplitude) can be immediately used 
for the next i + 1-th search. We shall show how this effect can be used to solve the problem of 
iterated search in a time O(^j^p) which is times as fast as by k sequential applications of 
Grover's algorithm. Thus we shall have a modification of the fast quantum search - the parallel 
quantum algorithm for iterated search which will be described later. In this article we take up 
mostly the particular case = 2 of IS, we name this problem repeated search (RS). 

RS-problem is connected with the known problem of structured search (SS). The problem of 
structured search is to find a unique solution xq, z/q of f{x, y) = I provided we have a function g 
whose support {x \ g{x) = 1} of cardinality M contains Xq. RS-problem is a particular case of 
SS when M = 1. The case 1 -C M <^ was investigated by Farhi and Gutmann in the work 
|[FG|| . They have found quantum algorithm for this case with a time complexity 0{\/MN), and 
also they wrote that the best known strategy for the case M = 1 is the sequential application of 
Grover algorithm. In the present paper it is shown how this evident strategy can be improved 
in constant factor a/2. At last note that our approach differs from the work |FG|. Farhi and 
Gutmann used only algebraic properties of Grover's algorithm whereas for RS-problem we need 
to work with an evolution of amplitudes in computation. 



3 Quantum computations and differential equations 



After early studies of R. Feynman (1^), P. Benioff ( [Pe|| ) and D. Deutsch ( 
approaches to quantum computations have appeared (look at |[BV|| , [ |Wa|| , ||Ya|| , 



numerous 
). Leaving 



aside the problem of decoherence and quantum codes (look at the articles [|AB|| , ||CLSZ|| , ||Ki2|| , 
pr| ) we shall regard ideal computations in closed systems. We use the abstract model of 
quantum computer independent of the formalism of classical algorithm theory. This model 
consists of two parts: classical and quantum (look at the picture 1). 
A state C of Classical part consists of the following objects. 

1) Registers with labels corresponding to transformations Ui- of the finite set {Ui} of ele- 
mentary unitary transformations with no more than 3 qubits each. (Strictly speaking transfor- 
mations on two qubits would suffice: look for example at the work |PiV|| ). Moreover, as follows 
from the works |n2|, | BBCDMSSSW ] there is a variety of possible choices of the set {Ui^}. 

2) Pointers aimed from these registers to as many qubits from the quantum part as there are 
arguments of the corresponding unitary transformation. Here each qubit is involved in exactly 
one transformation. 

3) Registers of an end of computation and of a query: e(C) and qu(C) respectively. 

A Quantum part is a tape partitioned into cells with one qubit each. Every qubit has two 
basic states |0), so its quantum state A|0) -|- |Ap -|- \fi\^ = 1, belongs to the curcle of 



radius 1 in two dimensional Hilbert space C^. 
belong to the tensor product H = (^C^ 



If n is a length of tape, all states of the tape 
...6^c' of spaces, corresponding to all qubits 



that is 7^ = C . Each state of quantum part is a superposition x 



N-l 

J2 AjCj of basic states 

i=0 
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N-1 

Co, ... ,6^-1 with complex amplitudes Aj where J2 l-^ip = 1, = 2"'. We can assume that all 

i=0 

eiG{0,ir. 

An Observation of this state x is a random variable which takes each value with the 
probability |Ajp. 

A Working transformation of quantum part corresponding to a fixed state of classical part 
has the form Uij^^Ui^^ . . Ui^ where each Ui. acts on qubits which the corresponding pointer 
aims to. 

Let /i, . . . , // be functions of the form {0, l}'^ — {0, l}"* and for each z = 1, 2, . . . , / there 
are special places in the quantum tape reserved for an argument of fi (query) and for a value 
of fi (answer). Denote by hi an initial contents of the place for answer. 

A Query transformation Quj is Quj^ (g) Quj^ (8> • • • (8> Qu/^ where for each i = 1, . . . , / 
Quj. Iflj, 6j) — > &i /j(aj)), is bitwise addition modulo 2. We name these functions Qu^, 
oracles. 

A Quantum algorithm is an algorithm determining evolution of the classical part: 

Co — ^ Ci — > . . . — > Ct (2) 

(in particular it determines a number T). A classical part plays the role of controller for 
quantum part and determines its evolution (look at the picture 2). 
A Quantum computation consists of two sequences: (2) and 

Qo Qi ^ • • • ^ Qt (3) 

where for each i = 0, 1, . . . , T — 1 e(Cj) = 0; eiCx) = 1 and every passage Qi — > Qi+i is: 

- a working transformation, corresponding to Ci, if qu(Cj) = 0, e(Cj) = 0, 

- a query transformation Quj, if qu(Ci) = 1, e{Ci) = 0. 

A result of this quantum computation is a contents of first qubits of quantum tape after 
the observation of final state Qt- An initial state (Co, Qo) of the computer is obtained from an 
input data x by some routine procedure. 

A time of computation (2), (3) is a number of query transformations (queries) in it. We see 
that in this model some oracles may be called simultaneously which causes interference between 
results of their actions. We shall prove that such interference can speed up computations in 
case of repeated search. 

From a physical standpoint systems of linear differential equations is a natural tool for the 
study of quantum computation. Wave function ip which determins an evolution of quantum 
computer satisfies Shrodinger equation ^ = iHip, where H is Hamiltonian within a real factor. 
An evolution of ip is determined by unitary operator U{t) : ip{t) = Uip{0) which satisfies the 
equation U = iHU . This equation is a prototype of all systems of differential equations studied 
below. We need only to choose Hamiltonian so that the amplitude of target state peaks in a 
point tquant which is less than the time of the best classical computation. In rare cases quantum 
parallelism makes it possible. 

Some other aspects of parallelism in computing may be found in ||BM|| , [PC)|| , ||LMS|| , |[Wo| . 
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4 An exact description of simple quantum search by 
differential equations 



4.1 Notations 

Assume the notations of Dirac where a vector a G as a column of coordinates is denoted 
by I a). A row obtained from \a) by the transposition and complex conjugation is denoted by 
(a|. A dot product of a,6 G C will be {a\b). A result of application of operator A to a vector 
\a) is denoted by A\a). For every transformations A, B of the forms Ci — > £2, ^2 — 
we denote by AB their composition which acts from right to left such that AB{x) = A{B{x)). 
Given vectors a ^ C, h E M. from linear spaces £, }A the state \a) <S)\h) E M. is denoted 
by |a,6). For a function F : F\X,Y) = \X,<p(Y)) we denote by F\y its restriction on 
Y: F\y\Y) = \ci>{Y)). 

Let / be a function of the form A — > A. We define an i-th iteration of /: /^*^ by the 
following induction on i. Basis: /^^^ = /. Step: /{'+^} = Z/^'^- 

4.2 Grover's quantum algorithm for simple search and its imple- 
mentation in our model 

Grover's algorithm for finding a unique solution xq of equation f{x) = 1 for a Boolean / is 
sequential applications of the following unitary transformation: G — —WRoWRt to the initial 

N-l 

state Xo = -j^ S where Walsh-Hadamard transformation \sW — J ^ . . J , 



|e), if 6 7^0, D u\ / if e 7^ Xq, 

~|0), if e = 0, "^'^^l \ -|xo), if e = xo. 

It is easily seen that W can be implemented on our model of quantum computer. 

N-l \n\-\i\ 

To implement Rt it is sufficient to apply Qu^ to the state ^ , where the 

last qubit is the place for oracle's answer. 

To implement Rq we need n + 1 ancillary qubits initialized by 0. Gonsider some function 
acting on three qubits: |main, ancilla, res) as follows. 



|000) - 


|000) 


|100) - 


— |ioi) 


1001) - 


1001) 


|101) - 


- Illl). 



Apply sequentially after each step moving pointers to right on one qubit in the main and 
ancillary areas (look at the picture 3). This makes res = 1 iff at least one of the main qubits 
is 1. Then inverse the phase of in the qubit res and fulfill all reverse transformations with 
in the reverse order restoring the initial states of ancillary qubits. 

Let Xi — (i'i |e') + &i|e), Xi+i — Gxi, e is a state of quantum part, corresponding to a 

target word Xq. The difference between xq and e is that e contains also ancillary qubits whose 
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values will be restored after each step of computation (it can be simply traced in what follows). 
The main property of Grover's transformation may be represented by the following equations 

(look at pBirq ] ). 




l-^)6. + 2(l-^ 



N 

-^bi + (1 - ^)ai 



(4) 



4.3 The passage to the system of differential equations 

The system (4) can be rewritten in the following form 




H, + 2(l-i)ai 



N 



(5) 



where the initial condition of quantum part gives ao = bo = l/yN. This system (5) with the 
initial condition is the system of difference equations approximating the following system of 
linear differential equations: 



-|6 + 2(l-^)a, 



(6) 



with two unknown functions a{t),b{t), constant 6 > and the initial condition a(0) = 6(0) = 
V^, where a^, 6, approximate a{i5), b{i6); 5 is a step. A difference between solutions of (5) and 
(6) on a segment of the form t G [0, 0{5\/N)] is 0{\/N 5^), hence the error of this approximation 
may be done as small as required by varying S (an integral part of number x is denoted by [x]). 
Solving (6) we obtain 



b+—-b + 0( 



0. 



6m 

Hence in within 0{^) a solution b of (6) can be approximated by a solution of equation 

4 



b + 



5^N 







(7) 



with the initial conditions 6(0) = |(6(0)5+ jjb(^)) = on the segment [0,2/ct;], where 

Lu = 2/6\/N. The required solution of (7) with this accuracy is 6 = sin(co't) + ;^cos(ti;t), and 

+ 

— 1 recurrent 



the maximum of amplitude 1 is in the point to = ^^-^5 — |. Then 

steps (4) are necessary and sufficient to achieve this value of b with this accuracy Thus we 









.<5. 
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obtain that the accuracy 0{-^) is reached in ^^^^ Grover's iterations. In the work |PBHT| 



its authors obtained an exact solution of (4) and thus proved that in fact a probability of 
error is even about The approximation of amplitude's evolution by systems of differential 
equations is more universal method. For example it makes possible to handle with the more 
involved case of parallel algorithm for RS which is the subject of next section. 
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5 Parallel algorithm for the repeated quantum search 



5.1 Definitions and result 



Let u, X, y he variables with values from three different copies of Tio = C 



N 



where Oi = 02 = :^(|0) — We assume the notations f2{x,y) for two oracles in the 

repeated quantum search and let ei, 62 be such values for x, y which represents unique solutions 
of equations /i = 1, /2 = 1. We denote the corresponding states of quantum tape by the same 
letters. 

Fut n = no^no^no^c\ Let 



Then 



Fi\u,x,y,a) 
F2\u,x,y,a) 
P\u,x,y,a) 



Fi\u,x,y,a) 



\u,x,y,ai® fi{u),a2), 
|M,x,?/,ai, 02 0/2(0;,?/)), 
|m0x, x,?/,a). 



\u,x,y,a), 
-\u,x,y,a), 



if u ^ Ci, 
if M = ei, 



F2\u,x,y,a) 



\u,x,y,a), 
-\u,x,y,a), 



if \x,y) 7^ |ei,e2), 
if \x,y) = |ei,e2); 



Define the following auxiliary unitary transformations on Ti: TZq = / Rqx Roy 1] VV = 

x,y,a2)P} where the lower indices x^y point the corre- 
sponding area of application for Walsh-Hadamard transformations and rotations of the phase 
of 0, J denotes the identity. 

The key unitary transformation of parallel algorithm for RS is 



(8) 



The parallel algorithm for RS is the sequential applications of Z beginning with the initial 

state 



N-l 



Xo = |0) 



EN) 



1 



1=0 



EN) 



i=0 



times. 



2V2 

Theorem 1 Let t = 

probability of error 0(-^). 



. Then the observation of Z^^^{xo) gives x = ei, y = 62 with 



5.2 An advantage of parallel quantum algorithm 

It follows from the definition of Z that oracles Fi and F2 for functions /i, /2 work simultaneously 



in parallel hence the algorithm requires approximately 



/N 



2V2 



simultaneous queries to obtain a 



result, when the sequential application of simple quantum searches with /i and then with /2 
time steps to obtain a result with the same probability. Note that for a simple 



requires 



search constant factor ^^^^ can not be essentially improved (look at |PBHT 



7 



Suppose that every query is fulfilled by a physical device (oracle) of the peculiar type 
corresponding to a form of query. A set with a minimum of oracles which is necessary for 
the solution of RS consists of one oracle for fi and one for /2. With these oracles we can 
run them simultaneously in the parallel algorithm and obtain a result y/2 times faster than 
by the sequential search. But if we have two copies of each oracle it is possible to achieve the 
same performance by sequential search if we divide the whole domain {0, 1}" into two equal 
parts of N/2 elements each and apply a simple quantum search at first with two copies of 
/i-oracle - one for each part, then, having ei, with two copies of /2-oracle. But this last way is 
expensive if every copy of oracle has a large cost, or impossible at all if every oracle is unique, 
say issues from a natural phenomenon. Just in this case of minimal possible set of oracles /i, 
/2 the application of parallel quantum algorithm for RS gives the increasing of performance 
in \/2 times. This speedup can be also obtained for IS if we apply this algorithm sequentially 
for the pairs i = 1,2,..., A; — 1. The resulting error probability will be 0{k/\/N). 

The remaining part of this work is devoted to the proof of Theorem and perspectives of this 
approach. 



5.3 A primary analysis of parallel algorithm for RS 

Note that each of Wy, Roy commutes with Wx, Rox, P, Fi, P commutes with F2, hence Z can 
be represented in the form 

Z = -[(/ (g) WxRoxWx (g) I)PF,P][-{I<^{WyRoyWy) (g) I)F2], 

or in the form 



Z = {-WxRo.WxJ',}{-WyRoyWyF2}, (9) 

where 

Ti\u,x,y,a) = 



\u, X, y, a) ii X ^ Ci, 
x, a) \i X = ei. 



The form (9) is exactly the repetition of Grover's transformations with oracles F2, J^i in this 
order, hence we can apply the formulas (4) for the resulting amplitudes of these transformations. 
Let the Z-iterations be Xo — > Xi — ^ ■ ■ ■ — ^ Xt, Xi+i = Z{xi), i = 0, 1, . . . , i - 1; 

Xi - h\eie2) + ai\eiN2) + ai\N^N2) + Al^^ie2), (10) 
where ei and 61,62 are the target states: unique solutions for fi{x) — 1 and for f2{x,y) — 1 

N 

respectively, A^i = e^, A^2 = J2 (we omit ancillary qubits). 

i=2 

We represent the transformation Xi — ^ ^iXi) ^ two sequential steps: 



X: 



i Zi{xi) ^x'i^ Mx'i) = Xi+i, 



where Zi = —WyRoyWyF2, Z2 — —WxRoxWxJ^i- To calculate the change of amplitude resulting 
from the application of Zi (or Z2) we shall fix a value of x (or y respectively). 

Step 1. Denote amplitudes of basic states in x'i by the corresponding letters with primes: 

X- = &:|eie2) + a',\e,N2) + a'^N,N2) + ^1^2). 
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Then for the two essentially different ways to fix a basic state for x: x = ei or x + Cj, j ^ 1 
we shall have the different expressions for new amplitudes. Use the property of the diffusion 
transformation WRqW to be an inversion about average (look at [ prl|| ). Let \av be an average 
amplitude of corresponding quantum state, 
a). X = Cj. 



A. 



iN-l)a,-bi 
N 



h' 



b). X 



X _ (Af-l)a,+ft , 



A' 



step 2. x'^ 

We have tw 
a), y 



_ 2(N~ 


-l)ai-2hi 




N 


_ 2{N- 


-l)ai-2bi 




N 


' 2Aa{; 




_ 2{Af 


-l)ai+2bi 




N 


_ 2(Af- 


-l)ai+2ft 




N 


) = Xi 


+1- 


ways 


to fix a 



- cti 
-(3i 



6,(l-|) + 2a,(l-^ 

-A, 

a,(l-f) + 2/34), 



-1 



62 or ?/ 



62- 



A -^"^ 


- m - b[ 


bi+i 


= 2A,, + b[ 




+ 2^(1 


N 




= 2Aa^, — A 


= m - It 


b). y = Cj, J ^ 2. 












A 






= 2Aai; + a'i 


= <(1 - If) 


+ 2a;.(l 


N 


"i+l 


2Aa-u ~ 




- ^) - a 



jV2. 



AT/' 



/_2_ 



AT/' 



Hence, the recurrent formulas for amplitudes of sequential steps 1,2 acquire the following 
form: 



= 6.(1 - f )2 + 2a,(l - i)(l - A) + 4a.(l - i)^ - 2^(1 - f )(1 - i); 
a,+i = a,(l - - 6,1(1 - I) + 2a,(l - |)(1 - i) + 2Af (1 - 

= «i(i - j^y + A^(i - ^) - - f )| + 6i]^; 
A+i = 2a.(l - i)(l - - A(l - 1)^ - k{l - |)f - 2a,(l - ^)|, . 

Thus the matrix of one step of algorithm has the form 

/I 2 4 -2 \ 
-^12^ 

V — N N 

A. _A 1 2l. 
\ ^ 2 -1 / 

The system of recurrent equations can be rewritten as the following system of difference 
equations. 

bi+i -h = 2a, + Aai - 2^ + 6,0i(^) + aiO^ij^) + a^O^ij^) + PiO^i^); 
a^+i-Oi = -j^bi + 2ai + ttiO^ij^) +biOG{^) + aiOrij^) + piOsi^); 
ai+i-ai = -j^ai + f^iOisi^) +aiOi4{^) + biOi5{^) +aiOiQ{^); 
A+i -Pi = -jfbi + 2ai - 2pi + a,Og{j^) + a,Oio(^) + AOii(^) + hOui^). 
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5.4 An approximation of amplitude's evolution by differential equa- 
tions 

Let {cj} be a sequence of vectors from C'^ : q = (cj, c?, . . . , cf), e C, which satisfies the 
following system of difference equations 

Cj+i - Cj = Aci, (12) 

where A is a matrix of size k x k with complex elements. 

Let m be an integer and a function C{t) : R — C'^ is a solution of the system of differential 
equations 

C{t) = mAC{t) (13) 

with the initial condition 

C(0) = Co. (14) 

Then the exact solution of the Cauchy problem (13), (14) is C{t) = R{t)co, where the resol- 
vent matrix R{t) = exp{mAt). The system (12) will be the system of difference equations 
approximating C{t) by Euler's method if we consider Cj as an approximation oi C{i/m), i — 
0, 1, . . . . The accuracy of approximation may be obtained by the Taylor formula C(^) = 
C(—) + —C(^) + TT-^Citi), < ti < Here the error ei of the one step of the rc- 

cursion (12) is the third summand 2^^(^'^) ~ ^A'^C{ti). Thus the error at the first step 
is ^A^ cxp{mA9i)co, at the second step: ^A"^ cxp{mA92)ci + exp{mA:^)^A^ ex]){mA9i)cO = 
^A^ exp{mA92){co + Acq + exp(mA6'i)co) + exp{A)^A'^ exp(mA9i)co, etc., at the i-th step 

the error will be Cj < | X) iexp(Aaj)A'^co, where < < 1. Hence if ||co|| < h, then the 

i=i 

error after i-th step is Cj = 0{ih). Particularly, for the initial conditions ||co|| = 0{j^) a good 
approximation can be obtained if i = o{N), and thus we can solve the Cauchy problem instead 
of (11) for I — 0{\/N) having error as small as required for sufficiently large A^. 

Define a new function B{t) as follows B{tm) — C{t). In terms of B the Cauchi problem 
(13), (14) acquires the form 

As(t)=AS(t), S(0) = Co. (15) 

Apply this to the solution q = a^, of (11), where cq = ;^)- Put B = 

\b,a,a,(3) for the scalar functions b,a,P,a and denote the argument of the function B by t. 
Then the equation (15) approximating (11) acquires the following form: 

b =2a + 4a-2/? + 60i(^) + ei + aOo(^); 

(16) 



a = 6 + 2a + e2 + 02(^)a 
P = -|6 + 2q;- 2/3 + 64 + 04(^)0; 



where Cj = aOoj(]^) + bOii{^) + P02i{j^) + aOsi{j^), i = 1,2, 3, 4, with the initial condition 

6(0) = a(0) = m = a{0) = 1. (17) 

Then for t = 0{VN), i = [t] the vector of error will be 5 = B{t) - Q = 0(1/ VN), N — > 00 
and with this accuracy we can write b{i) ~ bi for the amplitude bi of target state |ei, 62). 
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5.5 Tight analysis of the parallel quantum algorithm for RS 



Now we shall take up the system of linear differential equations (16) with the initial conditions 
(17). Our goal is to solve it on a segment of the form < t < 0{-\/N). The system (16) can be 
represented in the form B = MB, 
identity matrix) for the matrices 



where its matrix M = Z — l = Ari + E + H (1 denotes the 



/ 





2 


4 






2 

N 





2 










2 

N 








V 











J 



E 



f \ 





V 2 -2 y 



/ di di di 

(^2 d\ d\ 

(^2 '^2 d\ 

V d2 di di 



-2 + d^\ 
di 
di 
di J 



where di denotes different expressions of the form 0{N~ ), I = 1,2. 

We shall show that the deposit of Aq to the solution of (16), (17) is the main and deposits of 
E and H are negligible. What is the main difficulty here? Consider the resolvent matrix for the 
Cauchy problem (16), (17), this is the solution R{t) of the differential equation for matrices: 
R = MR with the initial condition -R(O) = 1. Then we have C{t) = RC{0). The matrix R 
has the form exp(Mt). But in our case the matrices Aq,E,H do not commutate, hence we 
cannot use the standard properties of exponent. In order to cope with this task we shall at 
first solve the Cauchy problem at hand neglecting deposits of E and H to the main matrix M. 
The legality of this approximation is shown in the Appendix. 

Now take up the reduced equation C{t) = AC{t) with the initial condition C(0) = cq. 
Excluding the last column and row containing only zeroes we obtain the new matrix Aq. The 



characteristic equation for is A + -|:A 



16 

Ar2 



and its nonzero solutions in within 0{j, 



N' 



are Ai_2 =~ 2\plij\fN . Then standard calculations give the approximation of solution as 



a 
a 



- 1 cos ^ 



;i8) 



J_ COS ^ + ^ 



irVN 



in within \0{^), 0{j^), 0{ ^^ )). The amplitude b from (18) peaks in the point ti — 
where b{ti) = 1 in within 0{-^). Assuming that deposits of E and H to the solution are small, 

steps of 



we obtain that the amplitude of target state |ei,e2) will be 1 — 0(l/vA^) after 
parallel algorithm which is \/2 times as small as the time of sequential quantum search. 



2^2 



5.6 Completion of the proof 

In the first version of this paper the deposit of E and H was estimated by the conventional 
procedure of approximation a solution of differential equation (look at the Appendix). This is 
the immediate but cumbersome way to prove that this deposit is vanishing. After the publi- 
cation of the first version of this paper Farhi and Gutmann informed me how to simplify this 
construction by uniting by pairs the sequential transformations in parallel algorithm (||FUT[). 
In this section I combine this idea with the approach of the first version. 
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At first turn to the orthonormal basis Ei = \e1e2), E2 = -^\eiN2), — j^\NiN2), 



E4 = ■;^\Nie2). The matrix Z in this basis acquires the form 

/ 



2^ 

N 



2 
N 



V 



N 
2 



4 

AT 

2 

2 1 _2_ 

1^-1} 

AT / 



Now group together each pair of unitary transformations in the algorithm: X2k — ^ X2k+i — ^ 
X2k+2 and denote by B the corresponding matrix: Bq — A\. It is sufficient to prove that 



r 7rv_iV 1 



|0,0, 1,0) — |1, 0,0, 0)11 = O(^), because one apphcation of Ai can only increase the 
error by O(^). 

The Cauchy problem for the recursion Cj+i = SoQ has the form c = {Bq — l)c, c = cq, and 
its resolvent has the form R — exp Bt, where B — Bq — 1. We have: 



B 



( 
-1 


V 



1 ^^ 
1 
-100 
/ 



in within 0{^). 

Thus we can consider only a projection of c to the subspace Ji\ spanned by Ei,E2,Es. 
Denote by D the matrix 








i 
V2 





\ 




i 





i 






V2 


V2 




v 





i 

V2 








Then the restriction of B to Hi has the form ^D. It is easily seen that D^''+'^ ^ D, D'^'' ^ 

then t = . Here in within we 



for k 
have 



1,2, 



If the number of steps is 

\b, a, a) 



2v/2 



exp(^Dt) = exp(7riL') = cos(7rD) + ism{7rD) 



1 - 



2 ' 4! 
1 -L)2(l -COSTT 



- . . . + i{nD - 



3! 



+ ... 



iL) sin TT = 1 - 2D^ 



that is 



/ 1 \ 

0-10. 

VI 0; 

The initial state is |0, 0, 1) in within Application of this matrix to the initial state gives 

|1, 0, 0) with this accuracy. Theorem is proved. 



6 Parallel implementation of iterated quantum search 
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6.1 Parallel quantum algorithm for IS 



Now take up an IS problem for arbitrary k. Consider an evolution of amplitudes arising 
when k oracles work in parallel. Let Xi = ^tol^i' ^2, • • • , N^) + a\\ei, N2, . . . , N^) + . . . + 
al\ei, 62, ... , Cfc) + Ri (it generalizes (10)), where Ri contains only basic states of the form 



.). The natural generalization of the transformation (8) will be 



W^i • • • /, each W-H transformation Wi acts on 
Rqi • • • (2) Rok 1, each rotation of O's phase i?oi 
• -^fc each Fi acts on Xi and inverses the sign of e. 



acts on a;j, t 



, I = 1,2, 
1, 2, . . . , fc, 
identities / act on 



where W^^^ ■■ 
. . . , k, TZq - 
= Fi 
ancilla. 

Let a matrix A determines an evolution of quantum state in parallel algorithm such that 
Xi = Axi-i- A represents the operator in 2'^"-dimensional space. We reduce A to an operator 
Ar acting on A; + 1-dimensional space generated by the vectors 

\Ni, N2, . . . , Nk), |ei, N2, . . . , Nk), . . . , |ei, 62, ... , Cfc). Then represent Ar as Ar = Aq + B, where 
Aq is Jacobi matrix of the form 



/ 



_2_ 

"TV 





V 



2 





2 







2_ 

' N 



\ 


2 
J 



(19) 



This matrix has the following nonzero elements: 2 - above the main diagonal and 
behind it. Its size is (fc + 1) x (A; + 1). Assume that the effect of reduction: A to Ar and the 
deposit of B are negligible. Generally speaking for k > 2 this assumption should be proved. An 
evolution of amplitude can be represented approximately as the solution of Cauchy problem 



Aa, a{0) = \N-''/\...,N-''/^), 



(20) 



where a{i) ~ \a\, . . . , a^), i integer, we assume that k <^ y/N. 

Given eigenvalues of matrix (19), a general solution can be obtained by the standard pro- 
cedure. These eigenvalues for Jacobi matrix is known (look at ||Bel|| , Chapter 2, ex. 32), there 
are Xm. = cos m9, m = 1,2, . . . ,k,9 = We shall not solve (20) here. 



6.2 Perspectives of parallel quantum algorithm for iterated search 

One can ask: can we obtain a speedup by a big constant factor for iterated quantum search 
when applying parallel action of more than two oracles? In all probability the answer is no. 

Estimate the growth of target amplitude from above. Canceling all — we can only increase 
this growth. This results in a simple system of linear differential equations whose solution is 

afc(t) = N-''/^ + 2^J N~^l^dt = N-^I^ + 2^N-^I\. The parallel algorithm for IQS can exceed 

sequential quantum search only if afc(t) is substantially large for t < k\/N. For example let the 
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parallel algorithm work only the quarter of time required for the sequential search. Then it can 

not reach the vanishing error probability for any fc, because for such t = '^^^ ~ ^ < 1. 

Nevertheless, parallel quantum algorithm has one advantage. Compare sequential and par- 
allel quantum algorithms for IS in case 1 -C A; -C \fN , if total time t = a/TV. If t = then 
t sequential applications of Z^^^ raise the amplitude up to the value ak(t) = a = 2 — e, 

where e — ^ (A^ — ^ oo). Hence the resulting probability is Ppar = • 

On the other hand if we have a total time then we can apply sequential quantum 
searches with the time \fN /k for each x^, i = 1,2, ... ,k. The probability Pgeg to find x° will 

be less than (^|^ because for one search it does not exceed . Consequently, Ppar exceeds 
Pseq in more than 2"^^ times. 

6.3 Conclusion 

To sum up, the parallel algorithm constructed above for repeated quantum search is a/2 times 
as fast as sequential application of fast quantum search and it requires the same hardware. 
The advance is taken of interference arising when two oracles act simultaneously on the set 
of entangled qubits. This parallel quantum algorithm can be applied to the problem of k 
dependent iterations of quantum search in areas of elements each, with the same effect of 
speedup in a/2 times. Here the error probability will be vanishing if A; = o{^/N), N — > oo. 

In addition, for the fixed total time ^/N a probability of success for parallel algorithm is 2"^^ 
times as big as for sequential algorithm. 

The effect of speedup in ^/2 times can not be increased essentially by the same procedure if 
we increase a number of oracles involved in the simultaneous action. Nevertheless, a possibility 
of further speedup of the iterated quantum search still remains. 
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A Appendix. An iterated approximation for solutions 
of differential equations 

A.l The method of simple iterations 

Here we present the first version of the proof. Let K, L, Kq be linear operators of the form 
C" — > C", K = Kq + L, a{t) : R — > C" be a vector function with scalar argument. 

Consider the following Cauchy problem for a system of linear differential equations with 
unknown vector function x{t) 

x{t) = Kx{t) + a{t), x(0) = xo e C™, (21) 

where Kq is regarded as the main summand and L determines a perturbation such that a 
solution Xo(t) of the problem 

X = KqX, x{0) = Xq 



17 



is known. Then a solution of (21) can be obtained as a limit of the sequence Xq, Xi, ... of vector 
functions depending on t, which are determined by the following recursion 

t 

Xi+i = xo + j {Kxi + a)dt. (22) 



Define the i-th difference between the approximations 1, 2, . . .. Then 

t t 
Si — J Lxodt, Si+i — J KSidt. Thus we conclude that 



t tl t2 ti~2 

Si^K' J J J ... J Si{ti_2)dti.2 . . . dt2dtidt. (23) 



A. 2 The method of complex iterations 

The simple method from the previous section can not be immediately applied to our situation 
because of inconvenient form of our main matrix M. We apply the following idea. In each 
step of approximation consider separately an approximations on the group of b, a, a and then 
on p. Here the main matrix A — Aq + B oi the system for b,a,a can be represented as a 
result of small perturbation B of almost nilpotent matrix which allows to estimate a deposit of 
perturbation. 

These matrices with the notations of (16) have the form 



/ 2 4 \ 

Ao= -I 2 

Consider the vectors 



/Oi(^) + On(^) Oo(^)+0oi(^) 03i(^) 

Ou{^) 02{^)+Oo2{^) 032(^) |. (24) 



= I - I + 014(^), 04(^) + Oo,{^), 2 + 034(^), -2 + 024(^)), 
It) =|-2 + 02l(i),022(^),023(i)). 

Preparing to define by induction a sequence of approximations to the solution of (16), (17) 
we introduce the notations 

Iq) \bii Cbi, Oii) , |Cj) \Ci, Pi) . 

Then the Cauchy problem (16), (17) is equal to the system 

\t) =A\c)+p\^), 

with the initial condition (17). 

Define the sequence of vector functions \ci) approximating the solution of (25), (17) by the 
following induction on i. 

Basis. |co) = \co,Po), where |co) = \bo,ao,ao) is a solution of Cauchy problem |c) = Ao|c), 
^(0) — If' ir) given by the approximations (18), Pq is a solution of equation Pq = (dplco) 
with the initial condition Po{0) = j^. 
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Step. Q+i is a solution of Cauchy problem 

Q+i = Aci+i + AI7/3), Q+i(0) = 1;^' ;^); 

is the solution of = (d/3|cj+i), /3j+i(0) = with the given vector function Cj+i obtained 
above. 

Put 5i = Ci — Ci-i = \S^, Sf, S", Si), i — 1,2, The main fact concerning these approxima- 
tions is the following 



Lemma 1 ForO<t< ^ 



Ni{i - 1)! ' ' - 7v^(i - 1)! ' "~ N-2+\i - 1)! 

This Lemma means that the vectors 5i form summable row with a module of sum less than 



20 

of matrix A. 



and thus we obtain Theorem. To prove the key Lemma 1 we need to clarify a main property 



A. 3 A nilpotency of the main matrix 

Denote by Cj any number whose absolute value is less or equal to (10/A^)-'. The main property 
of the matrix A is represented by the following Lemma. 

Lemma 2 For every j = 0, 1, . . . 



( ^i+i 

Cj+i Cj+i / V ^j+i 



^3 ^3 

A'^ = I 6,+i 6, 6, I , A'^^^ 



Proof 

Induction on j. 

A.4 Completion of the proof 

Proof of Lemma 1 

Induction on i. 
Basis. 

Consider the passage Cq — > Ci, consisting of two parts, we regard them as cases a) and b). 
a). Co — > ci. This passage is the sequence cq = — > c\ — > . . . — > Cq — > . . . — > ci, 

where c^"*"^ = c(0) + f{Ac?Q + f3o^)dt (7 in ordinary notations). The j-th difference will be 



A^. = cl- = / AA,^^dt, Ai = /(Ec^ + Pol)dt. 

Assume the following notation for an arbitrary function F{t): 
tW t t t _ _ 

J F{t)dt = J J ... J F{t)dt . . .dt dt. We will write a ^ 6 iff an absolute value of each 
000 
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component of a does not exceed the absolute value of the corresponding component of b. 
Then there exists such numbers Ig, s — 1, 2, 3 that for every k — 1,2,... we have / Ai :< 



t'' I h h h \ 

Now applying sequentially the equality (23) and Lemma 1, we obtain: 



r . , / lot \ 1 



Hence ci - Co = Ai + A2 + . . . ^ exp(^)|^, i, ^). 
b). ci — > ci. 

Introduce the following scalar function ^{b,a,a) = a, a). The difference for it is 

= - 

The first /3-difference 5f = /3i — /3o is the difference between the solutions of equations 
$0 = $(60, ao, ao) - 2/?o and /?! = $(61, oi, ai) - 2/3i. 

Lemma 3 For each i — 

t 

5f = exp(-2i) j ex.p{2t)A^idt. 


Proof 

Sf is the solution of the Cauchy problem + 25f = A$i, 5f (0) = 0, which can be obtained 
immediately. Lemma 3 is proved. 

Applying Lemma 3 we obtain |5f | < 2n-\/n ' Hence for the first step if i < y/N then 
Si Basis is complete. 

Step. 

The passage Cj_i — > Ci consists of two sequential steps: 

1) the passage from Ci_i to q, 

2) the passage from Pi — 1 to Pi. 

Take up Step 1. Applying the method of approximations (22) to the problem (25), we 
regard Cj as the limit of the sequence c^, c|, . . ., where = Q-i, 



t 







The accuracy of j-th approximation relative to the previous one is determined by the difference 
A^' = 1^^ ^j: ^j): given by 

A^- ^A^-'J A\t)dt, j>l, 

1 (26) 

A^ =JStildt. 
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Taking into account inductive hypothesis we obtain for the first difference 



^ / {i-iy.Ni/^ W-> w->w)^^- ^^^^ differences we have 



3*3- 



(i-l)W (i-l)!iV2' (i-l)W2 



Calculating the differences accordingly to the formula (26) and applying Lemma 2 we obtain 



+1 



(i-l)!Af'/2 
(i-l)!Ar'/2 



(2i-l)!Ar' (2i-l)!7V2 ' (2i-l)!Ar2 / ' 
SOt^J 10^ SOt^J 10^' 30t^^ 10^' 



(2j+l)!iVJ+i ' (2j+l)!Ari+2 ' (2j+l)!AfJ+2 



Summing over all j gives the inequalities for 5^, 5", 5° from the statement of Lemma 1. As for 
we estimate it as in the basis of induction, Then Lemma 3 gives \5f\ < jf^ij^]- Step of 
induction is complete. Lemma 1 is proved. Theorem is proved. 



21 



